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57] ABSTRACT 


Apparatus and methods are disclosed for providing fast 
decoding of Reed-Solomon and related codes. Cases of 
one and two data symbol errors are decoded directly 
from the remainder using a large pre-computed .table 
without calculating syndromes. Techniques for decod- 
ing cases of more than two errors are given where an 
optimized Chien search is used when more than four 
errors remain; when four or fewer errors remain, the 
Chien search is eliminated in favor of locating an error 
by direct solution of the error locator polynomial. The 
error locator and syndrome polynomials are adjusted 
after each error is found, and the error evaluator poly- 
nomial need not be computed. 


47 Claims, 21 Drawing Sheets 
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FAST REMAINDER DECODING FOR A 
REED-SOLOMON CODE 


BACKGROUND OF THE INVENTION 


This invention relates to information storage and 
retrieval systems, and more particularly to means for 
decoding codewords for use in error detection and 
correction in such systems. Even more particularly, this 
invention relates to Reed-Solomon codes and related 
codes, and to improved methods for the decoding of 
such codes. 

In a typical decoder, the coefficients of the syndrome 
polynomial, called frequency-domain syndromes, are 
first calculated by dividing the received codeword by 
each factor of the codé generator polynomial. This 
process is hardware-intensive and/or time-consuming, 
depending on implementation. 


It is possible to-efficiently generate the coefficients of ,, 


the remainder polynomial, called time-domain syn- 
dromes, by dividing the received codeword by the 
entire code generator polynomial. The time-domain 
syndromes contain the same information as the frequen- 
cy-domain syndromes. 

Methods of decoding using frequency-domain syn- 
dromes are known in the prior art, for example see 
Chapter 3 Practical Error Correction Design for Engi- 
neers by Neal Glover and Chapter 5 of Error-Correc- 
tion Coding for Digital Communications by Clark and 
Cain. For cases of one error, two frequency-domain 
syndromes can be used to generate an error location and 
value, which are then checked with other frequency- 
domain syndromes using a “Newton’s Identities” test. A 
similar procedure starting with four frequency-domain 
syndromes may be used to decode cases of two errors. 
The calculation of frequency-domain syndromes adds 
significant additional hardware and/or time delay to 
decoding of cases of one and two errors. 

For cases of more than two errors, the frequency- 
domain syndromes are typically used to iteratively gen- 
erate an error locator polynomial and an error evaluator 
polynomial. The error locator polynomial is typically 
used in a “Chien” search to locate the symbols in error, 
and the error evaluator polynomial is used to calculate 
the values of the symbol errors. This process is ineffi- 
cient when used to locate and evaluate errors in cases of 
less than five errors. 

It is thus apparent that there is a need in the art for 
improved decoding methods for Reed-Solomon codes 
and related codes. 


SUMMARY OF THE INVENTION 


Accordingly, it is an object of the present invention 
to provide improved decoding methods which decrease 
the correction time required for cases of one and two 
errors by decoding directly from time-domain syn- 
dromes without generating frequency-domain syn- 
dromes. ; : 

Another object is to provide improved methods to 
decrease the correction time required for cases of more 
than two errors by eliminating the Chien search when 
the number of yet-to-be-located errors is or has been 
reduced to four or less. 

Another object is to provide improved methods to 
decrease the correction time required for cases of more 
than two errors by calculating error values directly 
from the syndrome polynomial and error locator poly- 
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2 


nomial without generating the error evaluator polyno- 
mial. 

Another object is to provide improved methods to 
decrease the correction time required for cases of more 
than two errors by reducing the degree of the error 
locator polynomial and adjusting the syndrome polyno- 
mial as each successive error is located. 

Another object is to provide improved methods to 
decrease the correction time required for interleaved 
codewords of variable lengths by avoiding multiplica- 
tions by maintaining a parameter for computing point- 
ers to data symbols and remainder coefficients of the 
first such interleaved codeword and adjusting it as each 
successive codeword is corrected. 

The above and other objects are attained by the pres- 
ent invention which comprises means and methods for 
fast decoding of Reed-Solomon and related codes. 
Cases of one and two symbols in error are decoded 
directly from the read remainder, without computing 
syndromes. The location of a single symbol in error is 
quickly determined using the result of a single finite 
field division as an index into a small precomputed ta- 
ble. Using another larger precomputed table, the value 
is calculated and the validity of the location and value 
are confirmed. The larger precomputed table is also 
used to calculate and validate the locations and values 
of two symbols in error. Fast decoding of cases of more 
than two symbols in error is accomplished by employ- 
ing direct solution of an error locator polynomial -to 
locate errors when four or fewer errors remain to be 
located and by employing an optimized Chien search to 
locate errors when more than four errors remain to be 
located. The invention further simplifies and speeds the 
decoding of all cases of more than two symbols in error 
by eliminating the need to compute an error evaluator 
polynomial and by adjusting the error locator and syn- 
drome polynomials after each error is found. 

The above and other objects, features, and advan- 
tages of the instant invention will be more apparent 
from the following more particular description thereof 
present in conjunction with the accompanying draw- 
ings. 


BRIEF DESCRIPTION OF THE DRAWINGS 


FIG. 1 shows the environment in which the instant 
invention is used. 

FIG. 2 illustrates the steps required to decode Cases 
B, C.1, and C.2. : 

FIG. 2A illustrates the steps required to validate the 
error location and value for Cases B, C.1, and C.2 using 
sequential repeated blocks. 

FIG. 3 illustrates the steps required to calculate o; 
and o-2 for Case D using Ro through R3. 

FIG. 4 illustrates the steps required to decode Case A 
and calculate a and o2 for Case D using Rr.3 to Rr. 

FIG. 5 illustrates the steps required to calculate the 
error locations and values for Case D. 

FIG. 6 illustrates the steps required to validate the 
error locations and values for Case D. 

FIG. 6A illustrates the steps required to validate the 
error locations and values for Case D using sequential 
repeated. blocks. 

FIG. 7 illustrates the steps required to decode Case 
C.3. 

FIG. 8 illustrates the steps required to calculate the 
coefficients of S(x) required to decode cases of three 
and four errors for Case E. 
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FIG. 9 illustrates the steps required to iteratively 
generate the error locator polynomial for Case E. 

FIG. 9A illustrates the steps required to calculate the 
nth discrepancy d, using sequential repeated blocks. 

FIG. 10 illustrates the steps required to calculate the 
additional coefficients of S(x) required to decode cases 
of more than four errors for Case E. 

FIG. 11 illustrates the steps required to locate and 
evaluate errors by searching for roots of a(x) for Case 
E. 

FIG. 11A illustrates the modifications to FIG. 11 
required to evaluate a—/i*o-(x) at x=a/ using sequential 
repeated blocks. 

FIG. 12 illustrates the steps required to divide a(x) 
by (x@a“), compute the error value E, and adjust the 
coefficients of S(x) for Case E. 

FIG. 13 illustrates the steps required to transfer con- 
trol to the appropriate special error location subroutine 
for Case E. 

FIG. 14 illustrates the steps required to solve for a 
root X, and its log L, of a quadratic equation in a finite 
field. 

FIG. 15 illustrates the steps required to solve for a 
root X, and its log L, of a cubic equation in a finite field. 

FIG. 16 illustrates the steps required to solve for the 
log L of one of the four roots of a quartic equation in a 
finite field. 

FIG. 17 illustrates the steps required to correct errors 
in interleaved codewords. 


DESCRIPTION OF THE PREFERRED 
EMBODIMENT 


The following description is of the best presently 
contemplated mode of carrying out the instant inven- 
tion. This description is not to be taken in a limiting 
sense but is made merely for the purpose of describing 
the general principles of the invention. The scope of the 
invention should be determined with reference to the 
appended claims. 

In a finite field GF(2”), elements are composed of m 
binary bits and addition (@) consists of MODULO 2 
summation of corresponding bits; this is equivalent to 
performing the bit-wise EXCLUSIVE-OR sum of op- 
erands: 


x@y=xXORy. 


Note that subtraction is equivalent to addition since the 
MODULO 2 difference of bits is the same as their 
MODULO 2 sum. 

Multiplication (*) may be implemented using finite 
field logarithm and antilogarithm tables wherein LOG 
{a‘]=i and A LOG [i]=a!: 


x*y =0 
x*y = ALOG[LOG[x] + LOG[]] 


ifx =Oory=0 


ifx 4 O0andy +0 


where the addition of the finite field logarithms is per- 
formed MODULO 2”— 1. LOG [0] is undefined. 
Division (/) may be implemented similarly: 


x/y is undefined 
x/y =0 
x/y = ALOG[LOG[x] — LOGIy]] 


ify =0 
ifx = Oandy +0 


ifx ~Oand y+ 0. 
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4 
Note that for non-zero x, 
[x]=LOG [x] XOR 2”—1. 

In a decoder for an error detection and correction 
system using a Reed-Solomon or related code of dis- 
tance d for the detection and correction of a plurality of 
errors in codewords of n symbols comprised of 
n—(d—1) data symbols and d—1 check symbols, each 
symbol an element of GF(2”), a codeword C(x) is given 
by 


LOG [I/x]=—LOG 


C(x) =03—! Hx) @(x4—! #10) MOD G(X) (1) 


where I(x) is an information polynomial whose coeffici- 
ents are the n—(d— 1) data symbols and G(x) is the code 
generator polymomial 


ais Q) 
Ge) = 7 (K@aMOth © 


i=0 


where mois a parameter of the code. A code of distance 
d can be used to correct all cases of t= INT((d—1)/2) 
errors without pointers and is guaranteed to detect all 
cases of INT(d/2) errors. 

When e errors occur, the received codeword C’(x) 
consists of the EXCLUSIVE-OR sum of a transmitted 
codeword C(x) and the error polynomial E(x): 


C(x) =C(X)@E) (3) 
where 
E(x)=E, *xE1Q@ ... E,txLe (4) 


L; and E; are the locations and values, respectively, of 
the e symbol errors. 
The remainder 


R(x) =Rq—2*x4 2... BRi*X@RO (5) 
is given by 
R(x)=C(x) MOD G(x) (6) 


that is, the remainder generated by dividing the re- 
ceived codeword C’'(x) by the code generator polyno- 
mial G(x). 

By equation (1), 


C(x) MOD G(x)=0 : (7) 
so from equation (3), 
R(x)=E(x) MOD G(x) (8) 


For codes wherein d=6, the time required to decode 
cases of one and two errors may be reduced by decod- 
ing directly from the remainder by the use of a large 
pre-computed table fG,L). The value of each element of 
f(i,L) is the coefficient of the x! term of the remainder 
generated when x# is divided by G(x), that is, 


fid—2,L)*x4-2@ ... @f0,L)=x£ MOD G(x) (9) 


Note that fG,L)=40 for all i=0 to d—2 and L=d—1 to 
2m —2. 

Hereafter, expressions using elements of f(i,L) are 
written using the finite field antilogarithmic form of 
equation (9). Note however, that in the preferred em- 
bodiment of this invention the elements of f(i,L) are 
stored in finite field logarithmic form for computational 


5 


efficiency. Without loss of generality, discussion is lim- 


ited to codes wherein t2=3. To reduce the amount of 
storage required for the f(i,L) table, the one and two 
error correction procedure is structured to require only 
those elements of f(i,L) for L=d—1 to 2”—2 and i=0 
to T, where 


T=2t—1 for t=3, 
T=MAX(7,t+ 1) for t24, 


and only such elements are stored. 

From equations (4), (8), and (9), a single error at 
location L with value E will generate a remainder with 
coefficients 


Ri=EL) (10) 
when e errors occur in a codeword, the coefficients of 
the remainder are the EXCLUSIVE-OR sum of the 
coefficients of the remainders for each error taken sepa- 
rately: 


e qa)y 
R= 2 o BML). 
j= 


CASE A 


One or two check symbol errors 

A.1:,One error at location L<d—1 with value E. 
From equation (10), the remainder R(x) has a single 
non-zero coefficient 


Ri=E. 


A.2: Two errors at locations L}<d—1 and L2<d—1 
with values E; and E2. From equation (11), the remain- 
der R(x) has two non-zero coefficients 


Ri1E1, 
Rr2=EF2. 


' Cases A.1 and A.2 may be decoded by inspection of 
the remainder by counting the number of non-zero 
coefficients Rj. 

In the preferred embodiment of this invention, non- 
zero coefficients R; are counted for i=0 to t+1. If the 
total is less or equal to two, we are assured that there are 
at most d—t—1 errors, all at locations L;<d—1, none 
of which must be corrected since none affects a data 
symbol. The time required to decode cases of one or 
more data errors is reduced by counting non-zero coef- 
ficients only after a coefficient equal to zero has been 
detected in the course of decoding Cases B, C, or D. 


CASE B 


One error at location L=d—1 with value E. From 
equation (10), the coefficients of R(x) are given by 


Ri=EYGL). 


For given a and b, a=4b, the ratio 


flb,L) 
flaL) 


is unique for each value of L_ such that 
d—1SL<2"~—1. This yields a fast method for locating 
a single error. If either Rz=0 or Ry=0 there must be 
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6 

more than one error; otherwise the ratio R,/Ra may be 
used as an index into a pre-computed location table La,4 
derived according to 


BL 
is [2 | ats 


since from equation (10), 





Rp E*RbL) _“flbL) 
Ry E*ftaLl) ~ fal) ° 


Location table values corresponding to L<d—1 are set 
to an invalid number, for example zero, to flag the exis- 
tence of more than one error. 

In the preferred embodiment of this invention, the 
time required to determine the error location is reduced 
by using LOG [Rz/Rg] as an index into a location table 
Lg,p derived according to 


iz 


When the location L is known, the error value E may 
be trivially computed from equation (10): 


bL 


a [ rc [ 22 


Jal) 


Ra 
B= Fal) 
The computed error location and value may be rap- 
idly validated by gain using equation (10); test 


R;=E*fi,L) 


for values of i not equal to a or b; the indices of the 
coefficients R; used to determine the error location and 
value. i 

In the preferred embodiment of this invention, Ro and 
Ry are used to determine L and E, and coefficients R; 
are tested for i=2 to t. If all such coefficients are suc- 
cessfully tested, we are assured that there is one error at 
location LSd—1 with value E and at most d—t—2 
errors at locations t+1<Lj<d—1. 


CASE C 


Two errors, one error at location L;}2d—1 with 
value E; and a second error at location Lz<d—1 with 
value E2. From equation (11), the coefficients of R(x) 
are given by 


Ri=E 1 *KiL)), for iAL2; 


Ri=LE, FEL OLD, for i=L2. (12) 

In the preferred implementation of this invention, this 
case is divided into three subcases: 

C.1: t+1<L2<d—1 

Decoding proceeds as in Case B: error location and 
value computations using Ro and R; yield L; and E:, 
and coefficients R; are tested successfully for i=2 to t. 

C.2:2SL23t Error location and value computations 
proceed as in Case B, yielding L) and E}, but the test of 
Ry fails. If we perform one additional successful test of 
R:+1, we are assured that there is one error at location 
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L;Sd—1 with value E; and at most d—t—2 errors at 
locations 2<L;<d—1. 


o| = 


C.3: 0OSL2=1 

If either Ro or Rj is equal to zero, the procedure of 
Case B detects more than one error. If both Ro and Ri 
are non-zero, the error location and value computations 
of Case B yield incorrect results and more than one test 
of the coefficients R; fails. In either case we branch to 
Case D. 

If either Ro or Ri is equal to zero and the computa- 
tions of Case D are restricted to non-zero coefficients, 
the computation of the denominator of the coefficients 
of a(x) using Rr.3 to Rr will produce zero. If both Ro 
and R, are non-zero or the computations of Case D are 
not restricted to non-zero coefficients, the error loca- 
tion computations of Case D using Ro to R3 will cor- 
rectly locate the two errors, with one having location 
0=L;<1. In either case, we may then repeat the single 
error location and value computations of Case B but 
using R; and R;+1 in place of Ro and R, determining L 
using from a location table Ly:+; derived from 
f(t+1,L)/f(t,L)), and validating the results by testing 
the coefficients R; for i=0 to t—1, allowing one failure 
at i=0 or i=1. 

Case C.3 occurs with sufficiently low frequency that 
table storage space may be reduced without signficantly 
degrading performance by eliminating the L;;41 table 
and determining the location of the data symbol error 
by calculating L as a function of Ry41/Ri: 


(13) 
Rr+1 


G; 
a@ ae eae * R, 


tel 
G; 
Gry 


L=LOG 
Ri+t 


yom te 
Rr 





1g 
where G/' is the coefficient of the x‘ term of 


d-3 
G(x) = 7 (x @ a+) 
i=0 


In the preferred implementation of this invention, L 
for Case C.3 is computed from equation (13) using a=t 
and b=t+1. One validation testing failure is allowed 
for any i=0 to t—1. If validation testing succeeds, we 
are assured that there is one error at location L}>d—1 
with the value E;, a second error at location 
OSL2St—1, and at most d—t—3 errors at locations 
t+23L;<d-1. 


CASE D 


Two errors at locations Lj >d—1 and L22d—1 with 
values E; and E2. From equation (11), the coefficients of 
the remainder R(x) are given by 

RisE, PUL YBE2 ALL) (14) 

The coefficients 7; =a4!@a2 and o2=a£1*ql2 of 

the error locator polynomial o(x) may be computed 
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from four coefficients Rg, Rp, R,, and Rg according to 
the equations 


AabRaRb ® AacRaRe ® AadRaRa @ AbcRoRe @ AbdRoRa @ AcdReRd 


CabRaRs @ CacRaRe @ CadRaRd B CacRoRe B CadRoRa B CcaRcRa 


nan BabRaRs ® BacRaRe ® BadRaRd B BocRoRe @ BodRoRa O BeaRRa 
ae CabRaRy @ CacRaRe @ CadRaRa ® CocRoRe B CodRoRa @ CeaR-Ra 


20 


45 


50 


65 


where the pre-computed constants Ags, Bay, Cap, etc. 
are functions of a, b, c, and d of the form 


Aab = (a2? @ a24)*(G",_1*Ga" ® Ge"*G'"y_ 1), 
Bab = (a%*q25 ®@ at q2a#(G":_1*Gq" ® Ge"*G"y_}). 
Cab = (a7 @ a*)*(G",_1*Ga" @ Ge"*G"y_}). 


and G;’' is the coefficient of the x‘ term of 


d—4 ; 
G(x) = om (x @a™0t4, 
i=0 


with G_,” defined as zero. 

In the preferred embodiment of this invention, Ro to 
R3 are used in computing o-| and o2. For codes wherein 
t=4, the time required to compute a and ois reduced 
by choosing four non-zero zero coefficients of R(x): if 
any of Roto R3 is equal to zero, 0 and o-2 are computed 
using R7.3 to R7. If any of R73 to Rris also equal to 
zero, there exist more than two errors. 

We now have the coefficients of the error locator 
polynomial for the case of two errors 


x Qo} *x@o2=0. 


Solution of a quadratic equation in a finite field is 
known in the prior art; for example, see Chapter 3 of 
Practical Error Correction Design for Engineers by 
Neal Glover. Substituting x=y*o- yields 


a2 





xy Oy @c = 0, wherec = 


For each odd solution to this equation Y1, there is an 
even solution Y2=Y)@1. Y; can be fetched from a 
pre-computed quadratic table derived according to 


QUAD[? i] =i@ 1 for i=0, 2,... 2-2 


using c as an index. There are 2"—! such pairs of solu- 
tions; the other elements of the table are set to an invalid 
number, for example zero, to flag the existence of more 
than two errors. When Y;=40 and Y2=540 have been 
determined, reverse substitution yields expressions for 
the two error locations 


£1=LOG [o1*¥1] 
L2=LOG [o1 *¥2] 


If both L;<d—1 and neither L;is not equal to a or b, 
the indices of the coefficients R; used in Case B, there 
must exist more than two errors, otherwise the error 
situation of Case C.3 may exist. 

In the preferred implementation of this invention, 
Case C.3 is attempted if L;<d—1, whether or not L;is 
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equal to a or b; Case C.3 will detect cases of more than 
two errors. 
Solving the system of equations 


Ra=E, Kali @£2*flaL2) 


Ro= Ey FOL) BE2*KO.L2) 
yields expressions for the error values 


Ra*flb,L2) @ Ro*fla, £2) 
AGL AGL @ AGL AGL) 


Po Ra*flb,L1) @ Ro*haLi) 
2 RaLy ROL) @ AGL *Aa.L2) 

In the preferred embodiment of the invention, Ro and 
Rj are used in calculating error values if Ro to R3 were 
used in determining o(x), while R73 and Rr.2 are used 
in calcuating error values if R7.3 to Rr were used in 
determining o(x). 

The computed error locations and values may be 
rapidly validated using equation (11); test 


R=E PUL) BE2*fiL2) 


for values of i not equal to a, b, c, or d, the indices of the 
coefficients R; used to determine the error locations and 
values. 

In the preferred embodiment of this invention, coeffi- 
cients R;are tested for i=4 to t+ 1 if Roto R3 were used 
in determining locations and values, and coefficients R; 
are tested for i=0 to T-4 if R73 to Rr were used. If all 
such coefficients are successfully tested, we are assured 
that there are two errors at locations L;>d—1 and 
L2>d—1 with values E; and E2 and at most d—t—3 
errors at locations t+2SL;Sd—1. 


CASE E 


More than two errors, at locations L; with values E;. 

SYNDROME POLYNOMIAL GENERATION: 
We compute the coefficients S;of the frequency-domain 
. syndrome polynomial 


S(x)=Sq_2*x?-2@. .. BS *x@So 


from the coefficients R;of the remainder R(x) according 
to . 


d-2 
= = 


Si Rjtaimo+ i. 


Sequential computation of each coefficient S; would 
require d—1 references to each coefficient Rj. Physical 
constraints and interleaving of multiple codewords 
often make each reference to a coefficient Rj difficult 
and time-consuming. 

In the preferred embodiment of this invention, the 
time required to calculate the coefficients of S(x) is 
reduced by the following methods: 

(1) Each coefficient Rj; is referenced once and its 
contribution to each coefficient S; is sequentially 
computed and added. 

(2) A partial result for the contribution of each coeffi- 
cient Rj to each coefficient S;is maintained in finite 
field logarithmic form; 

(3) The use of software loops is minimized by using 
repeated sequential blocks. 
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(4) The calculation of coefficients S; required for the 
decoding of cases of more than four errors is de- 
ferred until the iterative error locator polynomial 
generation routine determines the existence of 
more than four errors. 

ERROR LOCATOR POLYNOMIAL GENERA- 
TION: The coefficients of S(x) are used to iteratively 
generate the coefficients of the error locator polynomial 
a(x). Such iterative algorithms are known in the prior 
art; for example, see Chapter 5 of Error-Correction 
Coding for Digital Communications by Clark and Cain. 

In the preferred embodiment of this invention, the 
time required to evaluate each successive discrepancy 
value is reduced by using the current degree of o-(x) as 
an index into a table of software jump addresses of 
appropriate starting points for each evaluation in a se- 
quence of repeated blocks. 

ERROR LOCATION AND EVALUATION: If the 
degree of o(x) indicates more than four errors exist, we 
evaluate o(x) at x=a# for each L, OSL <2”"~—1, until 
the result is zero, which signifies that a4 is a root of o-(x) 
and L is an error location. 

In the preferred embodiment of this invention, the 
time required to evaluate o(x) at each successive loca- 
tion is reduced by maintaining a software address 
pointer to the appropriate starting point for each evalu- 
ation in a sequence of repeated blocks, and by evaluat- 
ing wL*o(x) at x=a—/ where j is the degree of o(%), 
rather than evaluating o(x) at x=al. 

When the location L of an error has been abiennined, 
a(x) is divided by (x@a4), producing a new error loca- 
tor polynomial of degree one less than that of the old: 


— — 2%) 
ot) x@ak 4 

The error value E may be calculated directly from 
S(x) and the new o(x) using 


£ oP Sj-i 
E=aq—Lm0« i=0_ 
o(x) ab 
where j is the degree of the new a(x). 

In the preferred embodiment of this invention, the 
division of a(x) by (x@a£) and the calculation of the 
numerator and denominator of E are all performed in a 
single software loop. 

When the location L and value E of an error have 
been determined, the coefficients of S(x) are adjusted to 
remove its contribution according to 


S;=S;QE*a. L(m0+1) 


By reducing the degree of o(x) and adjusting S(x) as 
the location and value of each error are determined, the 
time required to locate and evaluate each successive 
error is reduced. 

When the degree of j of o-(x) is four or less, the time 
required to locate the remaining errors is reduced by 
using the special error locating routines below, each of 
which locates one of the remaining errors without using 
the Chien search. After the location of an error has been 
determined by one of the special error locating routines, 
its value is calculated o(x) is divided by (x@a“), and 
S(x) is adjusted in the same way as when an error is 
located by evaluating o(x). 
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When j=1, the error locator polynomial is 
x@o1=0 
By inspection, the root of this equation is @;=a“. Thus 5 
L=LOG [o}]. 


When j=2, the error locator polynomial is 





i 10 
x°Ba) *xGo2=0. 
One of the roots of this equation may be found using the 
quadratic method of Case D: 
15 
a2 
L=LoG [ erauap [=] 
2 
When j=3, the error locator polynomial is 20 
Poy *x?@o2*xOa3=0. 
Solution of a cubic equation in a finite field is known in 
the prior art; for example, see Flagg, U.S. Pat. No. ae 
4,099,162. Substituting 
x=w@o1, w=t@B/t, and v=t/B3/B 
yields a quadratic equation in v: 4a 
? @v@A3/B2=0 
where 
A=o}*@oz and B=o}*o2@o}. 35 
A root V of this equation may be found by the qua- 
dratic method above. Then by reverse substitution 
40 
= i A 
L=LOG [«: ® (BYV)3 ® Bn | 
When j=4, the error locator polynomial is 
45 
Qo} *PQo2*x?@o3*x@o4=0. 
Solution of a quartic equation in a finite field is known 
in the prior art; for example, see Deodhar, U.S. Pat. No. 
4,367,594. If o1=0, assign bj=07; for i=2 to 4, other- 50 
wise substitute 
a 
x® (03/01)! 
55 
to give 
A@br2z@b3zOHb4=0 
where 60 
by = i 
Fotos @ o1to2"73 @ oF 
53 = o1*b4, 65 


ba = ((oi*o3)! © 02)*b4. 


12 
The resulting affine polynomial may be solved in the 
following manner: 

GQ) Solve for a root Q of the equation 
g@b2*g@]b3=0 by the cubic method above. 

(2) Solve for a root S of the equation 
s?@b3/Q*s@b4=0 by the quadratic method above. 

(3) Solve for a root Z of the equation 22@O0*z@S=0 
by the quadratic method above. If 7; =0, L=LOG 
[Z], otherwise reverse substitution yields 


L=LOG [ even! @ +] 


Referring to FIG. 1, a data controller 100 having a 
host interface 102 is connected to a host computer 104. 
The data controller 100 also has a device interface 101 
which connects the data controller 100 to an informa- 
tion storage device 108. 

In the process of reading data from information stor- 
age device 108, data bits from information storage de- 
vice 108 are transferred through device information 
channel 116, through the device interface 101, and into 
the decoder 110. At the same time the data bits are being 
transferred into the encoder and time domain syndrome 
generator 110, they are transferred in parallel into the 
data buffer 106. As check bits are transferred into the 
decoder 110, remainder bits are generated by the de- 
coder 110 and transferred to the remainder buffer 107. 
After the data bits have been transferred into the data 
buffer 106 and remainder bits transferred to the remain- 
der buffer 107, the processor 105 uses the remainder bits 
from the remainder buffer 107 to detect and correct, if 
necessary, errors in the data bits in data buffer 106. 
After correction of any errors in the data buffer 106, the 
data bits are transferred through the host interface 102, 
through the information channel 118 through the host 
computer 104. 

Referring to FIG. 2, if either Ro or Rj is equal to zero, 
Step 200 transfers control to FIG. 4, having detected 
more than one error. Otherwise Step 202 determines the 
location £=£0,1{R1/Ro]. If L is invalid, that is, less than 
d—1, Step 204 transfers control to FIG. 3. Otherwise 
Step 206 calculates the error value E=Ro/f(O,L). Step 
210 prepares to validate L and E by initializing counters 
i=1, j=t—1, and k=2. Step 212 increments counter i. 
If R;-4E*f(i,L), Step 214 transfers control to Step 220. 
Otherwise Step 216 decrements counter j. If counter j is 
then: still greater than zero, step 218 transfers control 
back to Step 212 to continue testing coefficients Rj. 
Otherwise L and E have been validated and Step 250 
records them before transferring control to FIG. 17 for 
error correction. After a coefficient R; fails the test of 
Step 214, Step 220 decrements counter k. If counter k is 
then still greater than zero, Step 222 transfers control 
back to Step 212 to continue testing coefficients Rj. 
Otherwise L and E are invalid and control is transferred 
to FIG. 3. 

The time required for executing looping control and 
maintaining counters in FIG. 2 may be eliminated by 
replacing Steps 210 to 222 of FIG. 2 with the repeated 
sequential blocks of FIG. 2A. Steps 230 through 23x 
sequentially test coefficients R; for i=2 to t. If all of 
these tests succeed, L and E have been validated and 
Step 23x transfers control to Step 250 of FIG. 2. Other- 
wise control is transferred into Steps 240 to 24x, where 
sequential testing of coefficients R; continues through 


et Bee 13 


i=t+l. If all of the remaining tests in this second set . 


succeed, L and E have been validated and Step 24x 
transfers control to Step 250 of FIG. 2. If any of the 
tests of Steps 240 to 24x fails, L and E are invalid and 
control is transferred to FIG. 3. 

Referring to FIG. 3, on entry both Ro and Ri are 
known to be non-zero. If t24 and either R2 or R3 is 
equal to zero, Step 300 transfers control to FIG. 4 to 
attempt calculation of o and o2 using R7-.3 through Rr. 
Otherwise Step 302 calculates the products of Ro to R3 
taken two at a time. Step 304 computes the denominator 
of o-} and o>. If the denominator is equal to zero, Step 
306 transfers control to FIG. 8. Otherwise Step 308 
computes the log of the denominator, and the numera- 
tor of a}. If the numerator of 0°; is equal to zero, Step 
’ 310 transfers control to FIG. 8. Otherwise Step 312 
computes the log of o1, and the numerator of o>. If the 
numerator of o-2 is equal to zero, Step 314 transfers 
control to FIG. 8. Otherwise Step 316 calculates the log 
of o-2 and sets counter I=0 to signify that Roto R3 were 
used in computing o-1 and a2, then transfers control to 
FIG. 5. 

Referring to FIG. 4, on entry at least one of Ro to R3 
is known to be equal to zero. Step 400 counts the num- 
ber of non-zero coefficients R; for i=0 to t+1. If the 
result is less than three, there are no data symbols in 
error; Step 402 exits the correction procedure success- 
fully without correcting any data symbols, having de- 
coded Case A. Otherwise if any of R7.3 to Rrare equal 
to zero, Step 404 transfers control to FIG. 8, having 
detected more than two errors. Otherwise Step 406 
calculates the products of R73 to Rr taken two at a 
time. Step 408 computes the denominator of a and a2. 
If the denominator is equal to zero, Step 410 transfers 
control to FIG. 7, having possibly detected Case C.3. 
Otherwise Step 412 computes the log of the demonina- 
tor, and the numerator o}. If the numerator of o is 
equal to zero, Step 414 transfers control to FIG. 8. 
Otherwise Step 416 computes the log of o1, and the 
numerator of o-2. If the numerator of o-2 is equal to zero, 
Step 418 transfers control to FIG. 8. Otherwiste STep 
420 calculates the log of o2 and sets counter J=T—3 to 
signify that R7.3 to Rr were used in computing oj and 
o2, then transfers control to FIG. 5. 

In FIG. 5, which computes the error locations and 
values for Case D, Step 500 calculates the quantity 
c=o2/(o)2 and fetches ¥}=QUAD{c]; Y is known to 
be not equal to one, since o2 is known to be non-zero. If 
Y1 is equal to zero, Step 502 transfers control to FIG. 8. 
Otherwise Step 504 calculates locations L). If 1 is less 
than d—1, Step 506 transfers control to FIG. 7, having 
possibly detected Case C.3. Otherwise Step 508 calcu- 
lates Y2 and location L2; Y2 is known to be not equal to 
zero since Yj is not equal to one. If L2 is less than d—1, 
Step 510 transfers control to FIG. 7, having possibly 
detected Case C.3. Otherwise Step 514 computes the 
denominator of error values E) and Ep. If the denomina- 
tor is equal to zero, Step 516 transfers control to FIG. 8. 
Otherwise Step 518 calculates the numerator of E}. If 
the numerator of Ei is equal to zero, Step 520 transfers 
control to FIG. 8. Otherwise Step 522 calculates the 
numerator of Ep. If the numerator of E2 is equal to zero, 
Step 524 transfers control to FIG. 8. Otherwise Step 
526 calculates E; and E2, then transfers control to FIG. 
6. 

Referring to FIG. 6, Step 600 initializes counter i=4. 
If counter I is equal to zero, Step 602 transfers control 
to Step 606. Otherwise Step 604 initializes counter i=0. 
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Step 606 initializes counter j=t—2. If R=4E\*fEL1. 
J@L2*fG,L2), the computed error locations and values 
are invalid and Step 608 transfers control to FIG. 8. 
Otherwise Step 610 increments counter i and decre- 
ments counter j. If counter j is then still greater than 
zero, step 612 transfers control back to Step 608. Other- 
wise the error locations and values have been validated 
and Step 620 records them and transfers control to FIG. 
17 for error correction. 

The time required for executing looping control and 
maintaining counters in FIG. 6 may be eliminated by 
replacing FIG. 6 with thé repeated sequential blocks of 
FIG. 6A. Steps 650 through 65x sequentially test coeffi- 
cients R;for i=I to I1+-t—3. If any of these tests fails, the 
locations and values are invalid and control is trans- 
ferred to FIG. 8. Otherwise the locations and values 
bave been validated and Step 660 records them and 
transfers control to FIG. 17 for error correction. 

Referring to FIG. 7, if either R; or R;+1 is equal to 
zero, Step 700 transfers control to FIG. 8. Otherwise 
Step 702 calculates the location L from equation (13). If 
L is invalid, Step 704 transfers control to FIG. 8. Other- 
wise Step 706 calculates the error value E=R,/ftt,L) 
and prepares to validate L and E by initializing counters 
i=t, j=t—1, and k=2. Step 710 decrements counter i. 
If RiAE*fiL), Step 712 transfers control to Step 720. 
Otherwise Step 714 decrements counter j. If counter j is 
then still greater than zero, step 716 transfers control 
back to Step 710. Otherwise L and E have been vali- 
dated and Step 718 records them and transfers control 
to FIG. 17 for error correction. When a coefficient R; 
fails the test of Step 712, Step 720 decrements counter k. 
If counter k is then still greater than zero, Step 722 
transfers control back to Step 710. Otherwise L and E 
are invalid and control is transferred to FIG. 8. 

Referring to FIG. 8, Step 800 initializes all S;= Ro and 
initializes counters j=1 and JXM=mo. The purpose of 
counter JXM is to maintain the quantity j*mo for suc- 
cessive values of j by repaeated additions of mo and not 
by multiplication. If Rj=0, Step 810 transfers control to 
Step 850 since R, does not alter the coefficients S;. Oth- 
erwise- Step 820 computes k=LOG [Rj]+JXM and 
adds a* to So. Step 822 computes k=k-+ and adds a‘ to 
S1. Steps 824 to 8xx repeat the operations of Step 822 to 
compute and add the contribution of Rjto S;for i=2 to 
MIN(2t— 1,11). Step 850 increments counter j and com- 
putes JXM=JXM + mo. If counter j is then still less than 
d—1, Step 860 transfers control back to Step 810. Oth- 
erwise the coefficients S; for i=0 to MIN(2t—1,11) 
have been calculated and control is transferred to FIG. 
9. 

Referring to FIG. 9, Step 900 initializes the polyno- 
mials, parameters, and counters. Step 910 computes the 
nth discrepancy value dy. If dy, is equal to zero, Step 950 
transfers control to Step 960. Otherwise Step 952 up- 
dates o(x). If 1,21,+n—k, Step 954 transfers control 
to Step 958. Otherwise Step 956 updates o%(x) and 
other parameters. Step 958 updates o,(x). Step 960 
increments counter n. If n<t+1,, Step 970 transfers 
control to Step 972. If n5412, Step 972 transfers control 
back to Step 910 to perform another iteration; otherwise 
FIG. 10 computes the remaining coefficients S; before 
transferring control back to Step 910. Otherwise if 1, 
the degree of o(x), is greater than t, the number of 
errors the code can correct, Step 980 exits the correc- 
tion procedure unsuccessfully. If 1,=t, Step 980 trans- 
fers control to Step 982; if n=d—2, one additional itera- 
tion is required before terminating the algorithm so Step 
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982 transfers control back to Step 910. Otherwise we 
are assured that we have generated a valid error locator 
polynomial and control is transferred to FIG. 11. 

The evaluation of the nth discrepancy d, of step 920 
may be implemented with repeated sequential blocks as 
in FIG. 9A. Step 920 initializes d,=0. Step 930 uses 1, 
as an index into a table of software jump addresses to 
transfer control to the appopriate starting point in Steps 
940 to 94x for the evaluation of d,. Multiplying S, by 
0-70 is not required in Step 94x since oo is always equal 
to one. When d,, has been computed, control is trans- 
ferred to Step 950 of FIG. 9. 

Referring to FIG. 10, operation is similar to that of 
FIG. 8. After the remaining coefficients of S(x) have 
been calculated, control is transferred back to Step 910 
of FIG. 9. 

Referring to FIG. 11, Step 1100 initializes counters 
i=0, j=1,=the degree of o(x), and k=0=the number 
of error locations which have been found. If j is less 
than five, Step 1110 transfers control to FIG. 13. Other- 
wise Step 1130 evaluates a—/*o-(x) at x=al. If the 
result A is equal to zero, a root of o(x) has been found 
and Step 1160 transfers control to Step 1170, which 
forms L=—1i before transferring control to FIG. 12. 
Otherwise Step 1180 decrements counter i in the finite 
field. If counter i is then not equal to zero, Step 1190 
transfers control back to Step 1130 to evaluate o-(x) at 
the next location. Otherwise all possible locations have 
been tested without locating all the errors; therefore the 
correction correction procedure is exited unsuccess- 
fully. 

FIG. 11A shows modifications to FIG. 11 necessary 
for efficient evaluation of a—/*o(x) using repeated 
sequential blocks. Step 1120 is added to use counter j to 
compute an initial software jump address PTR=ADR(- 
Step 1140)+(¢—/)*K where K is the size of each Step 
1140 through 114x. Step 1130 initializes A=0, then Step 
1132 uses PTR to transfer control into Steps 1140 to 
114x, which sequentially evaluate a—/'*o-(x) at x=a! 
using Horner’s rule. Step 1170 is modified to update 
PTR=PTR+K to reflect the decrementing of j per- 
formed in FIG. 12. 

Referring to FIG. 12, Step 1200 increments counter 
k, records L, decrements counter j, then initializes D=1 
and N=§;. Step 1202 divides o(x) by (x@a4) and cal- 
culates the numerator and denominator of E£’- 
=E*ql*mQ, If the new ois equal to zero, the new o(x) 
has a root equal to zero, which is not the finite field 
antilogarithm of any error location; Step 1204 exits the 
correction procedure unsuccessfully. If the numerator 
is equal to zero, the computed error value is equal to 
zero, which is not a valid error value, so Step 1204 exits 
the correction procedure unsuccessfully. If the denomi- 
nator is equal to zero, the error value cannot be com- 
puted, since division by zero in a finite field is unde- 
fined; Step 1204 exits the correction procedure unsuc- 
cessfully. If oj, the numerator, and the denominator are 
each not eaqual to zero, Step 1204 transfers control to 
Step 1206, which calculates and _ records 
E=a-—LmQ*N/D. If counter j is equal to zero, Step 
1208 transfers control to FIG. 17 for error correction. 
Otherwise Step 1210 adjusts the coefficients of S(x) to 
remove the effects of the error just found and transfers 
control to FIG. 13. 

Referring to FIG. 13, if more than four errors remain, 
Step 1300 transfers control to Step 1180 of FIG. 11 to 
search for another error location by evaluating o-(x). If 
four errors remain, Step 1300 calls the quartic solution 
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subroutine of FIG. 16. If three errors remain, Step 1302 
transfers control to Step 1304, which sets parameters 
for and calls the cubic solution subroutine of FIG. 15. If 
two errors remain, Step 1302 transfers control to Step 
1306, which sets parameters for and calls the quadratic 
solution subroutine of FIG. 14. Otherwise one error 
remains and Step 1302 transfers control to Step 1308. 

If co is equal to zero, Step 1308 exits the correction 
procedure unsuccessfully, since the finite field loga- 
rithm of zero is undefined. Otherwise Step 1310 deter- 
mines L=LOG [o}] and transfers control to FIG. 12. 

Likewise, if one of the subroutines successfully deter- 
mines ‘an error location, Step 1320 transfers control to 
FIG. 12. Otherwise, the correction procedure is exited © 
unsuccessfully. . 

On entry to FIG. 14, the parameters c; and c2 de- 
scribe the quadratic equation 


x @Oey*¥xOeQ=0. - 


If c1=0, the equation has a repeated root. If c2=0, one 
of the roots is zero, whose log is undefined. If c1=0 or 
c2=0, Step 1400 exits the subroutine unsuccessfully. 
Otherwise Step 1402 determines a transformed root Y. 
If Y, is invalid, Step 1404 exits the subroutine unsuc- 
cessfully. Otherwise Step 1406 calculates the root X and 
its log L and returns successfully. 

On entry to FIG. 15, the parameters c}, c2, and c3 
describe the cubic equation 


3 @ey*x?@er*xOc3 =0. 


Step 1500 calculates the transform parameters A and B. 
If B is equal to zero, Step 1502 exits the subroutine 
unsuccessfully. Otherwise Step 1504 determines a root 
V of the quadratic equation 





Paya ae 
2 OBwW® am =0. 


using the QUAD table. If no such root exists, Step 1504 
produces zero and Step 1506 exits the subroutine unsuc- 
cessfully. Otherwise Step 1508 computes U. If U is not 
the cube of some finite field value T, Step 1510 exits the 
subroutine unsuccessfully. Otherwise Step 1512 calcu- 
lates T and a root X of the cubic equation. If X is equal 
to zero, Step 1514 exits the subroutine unsuccessfully. 
Otherwise Step 1516 calculates the log L of the root X 
and returns successfully. 

On entry to FIG. 16, the parameters 7}, 0-2, 073, and 
o-4 describe the quartic equation 


Bat 3Oo2*x2@o3*x@o4=0. 


If o1 is equal to zero, Step 1600 transfers contol to Step 
1610; if o3 is equal to zero, the quartic equation has 
repeated roots, so Step 1610 exits the subroutine unsuc- 
cessfully. Otherwise Step 1612 assigns b;=o7jfor i=2 to 
4 and transfers control to Step 1620. If oj is not equal to 
zero, Step 1600 transfers control to Step 1602, which 
calculates the numerator and denominator of transform 
parameter bq. If the denominator of b4 is equal to zero, 
Step 1604 exits the subroutine unsuccessfully. Other- 
wise Step 1606 calculates the transform parameters ba, 
b3, and b2 and transfers control to Step 1620. 

Step 1620 sets parameters for and calls the cubic 
solution subroutine of FIG. 15. If this returns unsuccess- 
fully, Step 1622 exits the subroutine unsuccessfully. 
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Otherwise Step 1630 assigns Q=X and sets parameters 

. for and calls the quadratic solution subroutine of FIG. 
14. If this returns unsuccessfully, Step 1632 exits the 
subroutine unsuccessfully. Otherwise Step 1640. sets 
parameters for and calls the quadratic solution subrou- 
tine of FIG. 14. If this returns unsuccessfully, Step 1642 
exits the subroutine unsuccessfully. Otherwise if o- is 
equal to zero, Step 1650 returns L successfully. Other- 
‘wise Step 1660 computes X. If X is equal to zero, Step 
1662 exits the subroutine unsuccessfully. Otherwise 
Step 1670 computes and returns L successfully. 

Referring to FIG. 17, on entry counter k is equal to 
the number of errors found, arrays L(*) and E(*) hold 
the pairs of error locations and values, D is the number 
of symbols in the data buffer and parameter DMN holds 
the reverse displacement plus one from the last symbol 
in the data buffer to the last data symbol in the current 
interleave. Step 1700 calculates the forward. displace- 
ment from the first symbol in the data buffer to the last 
data symbol of the current interleave. Step 1710 fetches 
the location L=L(k) of the next error to be corrected. 
If L is less than d—1, the symbol in error is a check 
symbol and need not be corrected, so Step 1720 trans- 
fers control to 1760. Otherwise Step 1730 calculates a 
forward displacement F within the data buffer from the 
reverse displacement L within the codeword interleave 
based on the data buffer length and number of inter- 
leaved codewords. If F is less than zero and therefore 
invalid, Step 1740 exits the correction procedure unsuc- 
cessfully. Otherwise Step 1750 corrects the data symbol 
in error with the error value E(k). Step 1760 decrements 
counter k. If counter k is then still greater than zero, 
Step 1770 transfers control back to Step 1710 to correct 
additional errors. Otherwise all errors have been cor- 
rected and the correction procedure is exited success- 
fully. 

While the invention herein disclosed -has been de- 
scribed by means of specific embodiments and applica- 
tions thereof, numerous modifications and variations 
could be made thereto by those skilled in the art with- 
out departing from the spirit and scope of the present 
invention. It is therefore to be understood that within 
the scope of the appended claims, the invention may be 
‘practised otherwise than as specifically described 

"herein. 
What is claimed is: : 
1. A decoder for an error detection and correction 
system using a Reed-Solomon code or related. code of 
degree d—1 for detection and correction of a plurality 
of errors in code words of n symbols comprised of k 
data symbols and d—1 check symbols, wherein each 
symbol is comprised of m binary bits of information and 
d, k, m, and n are positive integers, and further wherein 
t=INT((d—1)/2)23, said decoder comprising: 
data buffer means for storing said k data symbols; 
remainder generator means for dividing a codeword 
polynomial C(x) by a generator polynomial G(x) of 
said code and producing a remainder polynomial 
R(x) having remainder coefficients Ri; 

remainder buffer means for storing said remainder 
coefficients R; produced by said remainder genera- 
‘tor means; 

first table means comprising a table fG,L) wherein 

each element is comprised of error correction in- 
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18 
means for accessing said remainder buffer to re- 
trieve said remainder, 
means for applying said remainder to index said 
first table means and retrieve said correction 
information, and 
means for applying said correction information to 
_said data symbols in said data buffer to correct 
symbols that are in error. 

2. The decoder of claim 1 wherein each element of 
said first table means is comprised of a coefficient of the 
x‘ term of x MOD G(x). . 

3. The decoder of claim 1 wherein each element of 
said first table means is comprised of a finite field loga- 


- rithm of a coefficient of the x‘term of x MOD G(x). 


15 
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4. The decoder of claim 1 further comprising a sec- 
ond table means comprising a table Lg, wherein each 
non-zero element is defined as 


fb.L). 
Lab AGL) |- L 


wherein f(i,L) is a coefficient of the x/ term of x MOD 
G(x), L varies from d—1 to 2”"*—2 and further wherein 
O0<a<d—1, 0<b<d—l, and a+b. 

5. The decoder of claim 1 further comprising a third 


. table means comprising a table Lz,, wherein each non- 
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formation, and wherein OSi<d—1, and d—1<-_ 


L<2"—1; and 
processor means comprising 


zero element is defined as 


Lap] LOG [ 22- |]- L 


wherein f(i,L) is a coefficient of the x/term of x MOD 
G(x), LOG [f(b,L)/f(a,L)] is the finite field logarithm of 
f(b,L)/f(a,L), L varies from d—1 to 2”—2 and further 
wherein 0<a<d—1, 0Sb<d—1, and a+b. 
6. The decoder of claim 5 wherein said means for 
applying said remainder comprises: 
means for determining a location L of a data symbol 
error; 
means for calculating a value E of said data symbol 
error; and 
means for validating said location L and said value E 
of said data symbol error. 
7. The decoder of claim 6 wherein said means for 
determining said location L comprises: 
means for calculating an index value 


Ro 
LOG zz f 


wherein Rg and R,z are two of said coefficients Rj 
means for applying said index value to said third table 
means to retrieve said location L of said data sym- 

bol error. 
8. The decoder of claim 6 wherein said means for © 
determining said location L comprises: 
means for calculating L according to 


(13) 
G'a 
G'b 


Re 
Ra 


* 








a? @al* 
L=LOG 
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wherein G/ is a coefficient of the x/ term of 


d-3 
G(x) = 7 (x @a+h, 


i= 


9. The decoder of claim 6 wherein said means for 
calculating said value E comprises: 
means for applying said location L to produce an 


index value for element f(a,L) of said first table 10 


means; 
means for applying said index value to said first table 
means to retrieve said element f(a,L); and 
means for dividing a coefficient Rg by said element 
f(a,L) to produce said value E of said data symbol 
error. 
10. The decoder of claim 6 wherein said means for 
validating said location L and said value E comprises 
means for calculating a finite field logarithm of said 
value E; . 
means for testing a plurality t of said remainder coef- 
ficients R; wherein isa and i=+b, each test compris- 
ing a sequential repeated block comprising: 
means for adding said finite field logarithm of said 
error value E to an element f(i,L) of said first 
table means to produce a finite field logarithm of 
a test value; 
means for calculating a finite field antilogarithm of 
said finite logarithm of said test value; and 
means for comparing said finite field antilogarithm 
_ of said test value to said coefficient Rs 
means for counting a number of said tests wherein 
said finite finite field antilogarithm of said test 
value is not equal to said remainder coefficient Ry 
means for recording an indicium when said number is 
less than two; and 
means responsive to said indicium for using said loca- 
tion L and value E to correct said data symbol 
error. 
11. The decoder of claim 5 wherein said means for 
applying said remainder comprises: 
means for computing parameters oj and 2; 
means for determining locations L; and L2 of two 
data symbol errors from oj and 07; 
means for calculating values E; and E2 of said two 
data symbol errors; and 
means for validating said locations and said values of 
said two data symbol errors. 
12. The decoder of claim 11 wherein said means for 
computing said parameters o- and o2 comprises: 
means computing non-zero parameters D, Ni, and 
N2 according to 


Ny = AabRaRb ® AacRaRe ® AadRaRa @ AacRoRe @ AsaRoRa @ 
AcdR-Ra 

N2 = BabRaRs ® BacRaRe ® BadRaRd ® BocRoRe ® BodRoRa ® 
BedReRa 

D = CapRaRp ® CacRaRe ® CadRaRa @ CocRoRe ® CodRsRa O 
CedReRd 
wherein pre-computed constants Ags, Aac, Aad, 


Abes Abd; Acdy Bab, Bae, Baa; Boe Bod, Bea Cab, Cac, 
Cad, Coc, Cog, and C,g are functions of a, b, c, and d 
given by: 
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dab = (024 @ a*)*(G"e_\*G"d ® G"e*G"d—1), 

Age = (074 @ a*)*(G"p_1*G"a ® G"s*G"d—1), 

Aad = (a4 @ @4*(G"5_1*G"e B G'5*G'e— 1), 

Abe = (a? & a)*(G"a—1*G"d © G'a*G"g—1). 

Abd = (a? @ a?4)*(G"q_1*G"e @ G'a*G'e—1)s 

Acd = (a @ 0*4)*(G"g_1*G"b © G'g*G"5—1). 

Bap = (ata? @ abta?)*(G"c_\*G"d ® G'*G"'a—1) 
Bac = (a%*a* © a°ta”)*(G"5_1*G"a © G"y*G"ad— 0), 
Bad = (a%*a?4 @ af*a29)*(G"4_|*G"e @ G'5*G"e_1), 
Bre = (ata @ a°ta?)*(G"a_1*G"'d ® G'a*G"d—1)s 
Boa = (abta?4 @ at*a?)*(G"g_1*G"c ® G'a"G'e— 1) 
Bed = (a°ta*4 ® a#a)*(G"g_1*G"s ® G'a*G"b—1) 
Cap = (a7 @ a*)*(G"e-1*G"a ® G'FG"a_1), 

Cac = (0? © a)*(G"b—1*G'd ® G"s*G"'a—1), 

Cad = (a? © a4)*(G"5_1*G"c ® G'5*G"e—1), 

Cte = (ab ® a)*(G"a_1*G'a ® G'a*G"a—1), 

Cod = (a° ® a4)*(G"a—1*G"e © G'a*G"e_1), 


Cea = (0° @ a4)*(G"g_1*G"b © G"a*G"5-1), 


and G;’ is a coefficient of the x‘ term of 
d—4 


am (x @amt+h 


i= 


G"(x) = 


with G_;” defined as zero; and 

means for computing o1;=N;/D and 02=N2/D. 

13. The decoder of claim 11 wherein t24 and se- 
lected coefficients Ra, Rp, Re, and Rgof said coefficients 
R; are each not equal to zero. 

14. The decoder of claim 11 wherein said means for 
calculating said values E; and E2 of said two data sym- 
bol errors comprises: 

means for computing non-zero parameters D, Ni, and 

N2 according to 


D = flaLi)*fb.L2) ® fbL)*fla£2); 
M = Ra*fb.L2) ® Ro*fha.L2), 


N2 = Ra*fl,L1) ® Re*fla.L1); and 


means for computing E;=Nj/D and Ez=N?/D. 
15. The decoder claim 11 wherein said means for 
validating said locations and said values comprises: 
means for calculating finite field logarithms of said 
values E; and E2; 
means for testing a plurality t—2 of said remainder 
coefficients R; wherein is4a, is¢b, is4c, and is4d, 
each test comprising a sequential repeated block 
comprising: 
means for adding said finite field logarithms of said 
values E; and E2 to respective elements f(i,L1) 
and f(i,L2) of said first table means to produce 
finite field logarithms of test values, 
means for calculating finite field antilogarithms of 
said finite field logarithms of said test values; 
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means for calculating an EXCLUSIVE-OR sum of 
said finite field antilogarithms of said test values; 
and 
means for comparing said EXCLUSIVE-OR sum 
of said finite field antilogarithms of said test val- 
ues to said coefficient Rs 
means for recording an indicium when each said 
EXCLUSIVE-OR sum of said finite field antiloga- 
rithms of said test values is equal to said coefficient 
Rg and : 
means responsive to said indicium for using said loca- 
tions and values to correct said data symbol errors. 
16. The decoder of claim 1 wherein said means for 
applying said remainder comprises: 
means for counting a number of non-zero coefficients 
Ryin a plurality of t+ 1 of said coefficients R;in said 
remainder buffer and recording an indicium when 
said number is less than three; and 
means responsive to said indicium for terminating 
error correction successfully. 
17. The decoder of claim 1 wherein said means for 
applying said remainder comprises: 
means for validating locations Lj and values E; of e 
errors comprising means for testing a plurality of 
said remainder coefficients R; according to the 
equations 


Ri= ee o Et fi.L). 

18. The decoder of claim 1 wherein said processor 

means comprises: 

means for computing a syndrome polynomial S(x) 
from said remainder polynomial R(x); 

means for generating an error location polynomial 
o(x) from said syndrome polynomial S(x); 

means responsive to said error locator polynomial 
o(x) for locating errors; and 

means responsive to said error locator polynomial 
o(x) and said syndrome polynomial S(x) for evalu- 
ating errors. 

19. The decoder of claim 18 wherein said means for 

computing said syndrome polynomial comprises: 

(1) means for initializing a coefficient So and all other 
coefficients S; of said syndrome polynomial S(x) to 
a coefficient Ro of said remainder polynomial R(x); 

(2) means for initializing a counter j to 1; 

(3) means for computing a finite field logarithmic 
partial result comprising a MODULO 2”—1 sum 
of a finite field logarithm of a non-zero coefficient 
Rj and j*mo; 

(4) means for calculating a finite field antilogarithm of 
said partial result; 

(5) means for EXCLUSIVE-OR adding said finite 
field antilogarithm of said partial result to said 
coefficient So; 

(6) means for MODULO 2”~— | adding said counter j 
to said partial result; 

(7) means for calculating a finite field antilogarithm of 
said partial result; 

(9) means for EXCLUSIVE-OR adding said finite 
field antilogarithm of said MODULO 2”—1 sum 
and one of said coefficients S; 


(10) means for repeating said means (6) through (9) 
for said coefficients S; wherein i=2 to 
MIN(d—2,11); 


(11) means for incrementing said counter j; and 
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(12) means for repeating said means (3) through (11) 
for said coefficients Rj wherein j=2 to d—2. 

20. The decoder of claim 18 wherein said means for 

generating said error locator polynomial comprises: 

means for computing an nth discrepancy d, compris- 
ing sequential repeated blocks, each said block 
comprising means for calculating a MODULO 
2™—1 sum of said nth discrepancy d,, and a finite 
field product ‘of a coefficient op; and a coefficient 
Sn-6 

means using a degree |, of the nth error locator poly- 
nomial as an index into a table of software ad- 
dresses of each said block; 

means for recording an indictum when said parameter 
n is equal to twelve; and 

means responsive to said indicium for generating 
coefficients S; for i=11 to d—2. 

21. The decoder of claim 18 wherein said means for 

locating said errors comprises: 

means for recording an indicium when said error 
locator polynomial o(x) is of degree j greater than 
four; 

means responsive to said indicium for locating one of 
said errors comprising evaluating said error locator 
polynomial o(x) for successive values of L until 


A = dlk*g(x)|qg—L = 0, 


said means for evaluating comprising a sequence of 

repeated blocks each said block comprising means 

for calculating a finite field product of a—4 and an 

EXCLUSIVE-OR sum of said parameter A and a 

coefficient oof said error locator polynomial o-(x); 

and 

means for maintaining a software address of a starting 
point for next said evaluation. 

22. The decoder of claim 18 wherein said means for 

locating said errors comprises: 

means for recording an indicium when said error 
locator polynomial o-(x) is of degree j less than or 
equal to four; and 

means responsive to said indicium for locating one of 
said errors comprising: 

means for calculating a finite field logarithm of a 
root of a quartic equation in a finite field; 

means for calculating a root and a finite field loga- 
rithm of said root of a cubic equation in a finite 
field; 

means for calculating a root and a finite field loga- 
rithm of said root of a quadratic equation in a 
finite field; and 

means for calculating a finite field logarithm of a 
root of a linear equation in a finite field. 

23. The decoder of claim 18 wherein said means for 

evaluating said errors comprises: 

means for dividing said error locator polynomial o-(x) 
by (x @a#) to produce a new error locator polyno- 
mial o(x) and calculating an error value E from 
said syndrome polynomial S(x) and said new error 
locator polynomial o(x), comprising a single soft- 
ware loop comprising: 

(1) means for initializing a counter g=1, a remain- 
der R=1, a denominator D=1, and a numerator 
N=o0; 

(2) means for calculating a MODULO 27'—1 sum 
of said remainder R and a finite field product of 
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a finite field antilogarithm of said location L and 
said remainder R; 

(3) means for storing said MODULO 2”'—1 sum as 
said remainder R and as a coefficient og of said 
error locator polynomial o(x); 5 

(4) means for calculating a MODULO 2”—1 sum 
of said remainder R and a finite field product of 
a finite field antilogarithm of said location L and 
said denominator D; 

(5) means for storing said MODULO 2”— 1 sum as 
said denominator D; 

(6) means for calculating a MODULO 2”—1 sum 
of said numerator N and a finite field product of 
said remainder R and a coefficient Sj, of said 
syndrome polynomial S(x); 

(7) means for storing said MODULO 2”— 1 sum as 
said numerator N; 

(8) means for incrementing said counter g; and 

(9) means for repeating said means (2) through (8) 
for values of said counter g up to and including 
J 

means for recording an indicium when R, D, or N is 
equal to zero after the operation of said means (1) 
through (9); 

means responsive to said indicium for terminating 
error correction unsuccessfully; 

means for calculating a finite field quotient of said 
numerator N and said denominator D; 

means for recording a finite field logarithm of said 
finite field quotient as a parameter E’; 

means for calculating a finite field product of said 
finite field quotient and a finite field antilogarithm 
of —L*mg; 

means for recording said finite field product as said 
error value E; and 

means for adjusting coefficients of said syndrome 
polynomial S(x) comprising a software loop com- 
prising: 

(a) means for initializing a counter g=0; 

(b) means for calculating a finite field antilogarithm 
of said parameter E’; 

(c) means for calculating a MODULO 2”—1 sum 
of said finite field antilogarithm and a coefficient 
S; of said syndrome polynomial S(x); 

(d) means for storing said MODULO 2” — 1 sum as 
said coefficient S; 

(e) means for calculating a MODULO 2”—1 sum 
of said parameter E’ and said location L; 

(f) means for storing said MODULO 2” — 1 sum as 
said parameter E’; 

(g) means for incrementing said counter g; and 

(h) means for repeating said means (b) through (g) 
for values of said counter g up to and including 
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24. The decoder of claim 1 wherein m=8, t=8, G(x) 55 
is a GF(256) polynomial 


15 ; 
Gx) = 7x Ga, 
0 


es 60 
and a/are elements of a finite field generated by a GF(2) 
polynomial 

8Qxt@Oxr@xOl. 65 


25. The decoder of claim 1 wherein m=8, t=8, G(x) 
is a GF(256) polynomial 


24 


15 
Gx) = 7 (x @ a4, 
i=0 


a 
mo= 120, and a’ are given by 
a! =(beta‘)88, 


wherein beta! are elements of a finite field generated 
* by a GF(2) polynomial 


BOs @rOr@l. 


26. The decoder of claim 1 wherein m=8, t=4, G(x) 
is a GF(256) polynomial 


7 
Gx) = 7e@ a0+ fy, 


i= 


mo= 124, and a/ are given by 
a!=(beta‘)88, 
wherein beta! are elements of a finite field gener- 
ated by a GF(2) polynomial 


BQ OrgrOl. 


27. The decoder of claim 1 wherein m=8, t=2, G(x) 
is a GF(256) polynomial 


3 
Gx) = 7 (x @ ath, 


i= 
mo= 126, and a/are given by 
a! =(betal)88, 


wherein beta! are elements of a finite field gener- 
ated by a GF(2) polynomial 


BQ i Qv@xrol. 


28. In a decoder for an error detection and correction 
system using a Reed-Solomon code or related code of 
degree d—1 for detection and correction of a plurality 
of errors in codewords of n symbols comprised of k data 
symbols and d—1 check symbols, wherein each symbol 
is comprised of m binary bits of information and d, k, m, 
and n are positive integers, and further wherein 
t=INT((d—1)/2)23, an error decoding method com- 
prising the steps of: 

storing said k data symbols in a data buffer; 

generating a remainder polynomial R(x) having re- 

mainder coefficients R; by dividing a codeword 
polynomial C(x) by a generator polynomial G(x) of 
said code; 

storing said remainder coefficients in a remainder 

buffer; 

applying said remainder coefficients to index a first 

table fG,L), each element of said table being com- 
prised of a coefficient of the x/ term of x MOD 
G(x) wherein L varies from d—1 to 2’"*—2 and i 
variés from 0 to d—2, to produce correction infor- 
mation; 6 
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applying said correction information to said data 
symbols in said data buffer to correct symbols that 
are in error. 

29. The method of claim 28 wherein said step of ap- 
plying said remainder coefficients further comprises the 5 
steps of: 

counting a number of non-zero coefficients R; in a 

plurality t+2 of said coefficients R; in said remain- 
der buffer; and 

terminating error correction successfully when said 

number is less than three. ; 

30. The method of claim 28 wherein said step of ap- 
plying said remainder. coefficients further comprises the _ 
steps of: -_ 

determining a location L of a data symbol error; 

calculating a value E of said data symbol error; and 

validating said location L and said value E of said 
data symbol error. 

31. The method of claim 30 wherein said step of de- 
termining said location L comprises the steps of: 

‘calculating an index value from two coefficients Rg 
and Ry, of said coefficients R;: 


Rp 
LOG Ri and 


applying said index value to a second table La» 
wherein each non-zero element is given by 


Eg [ +03 [#2 |]- Zz 


for L=d—1 to 2™—1 to retrieve said location L of 35 
said data symbol error. 
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dividing a coefficient Rg by said element f(a,L) to 

produce said value E of said data symbol error. 
34. The method of claim 30 wherein said step of vali- 
dating said location L and said value E comprises the 
steps of: 
calculating a finite field logarithm of said value E; 
testing a plurality t of said remainder coefficients R; 
wherein is4a and i=4b, each test comprising a se- 
quential repeated block comprising the steps of: 
adding said finite field logarithm of said error value 
E.to an element f(i,L) of said first table to pro- 
duce a finite field logarithm of a test value; 
calculating a finite field antilogarithm of said finite 
field logarithm of said test value; and 
comparing said finite field antilogarithm of said test 
value to said coefficient R; ; 
counting a number of said tests wherein said finite 
finite field antilogarithm of said test value is not 
equal to said remainder coefficient Rj; and 

correcting said data symbol error using said location 
L and said value E when said number is less than 
two. 

35. The method of claim 28 wherein said step of ap- 
plying said remainder coefficients further comprises the 
steps of: : 

computing parameters oy and o>; 

determining locations L; and Lz of two data symbol 

errors from o-; and 07; 

calculating values E; and E2 of said two data symbol 

errors; and 

validating said locations and said values of said two 

data symbol errors. 

36. The method of claim 35 wherein said step of com- 
puting said parameters o) and o-2 comprises the steps of: 

computing non-zero parameters D, Nj, and N2 ac- 

cording to 


Ni = AabRaRo @ AacRaRe @ AadRaRd B AvcRoRe BO AsdRoRa B AcadRcRa 


N2 = BabRaRo @ BacRaRe @ BadRaRa B BocRsRe B BodRoRa OB BedReRa 


D=CaoRaRb B CacRaRe @ CadRaRa DB CocRoRe B CodRoRa OB CedRcRa 


32. The method of claim 30 wherein said step of de- * 
termining said location L comprises calculating L ac- 
cording to 





. (13) 
a bs Ga * Ro a 
at @a Ge Ra 
L= LOG ; 
i Ga, Re 
O-Gy "Ra 
. . . ry 55 
wherein G;’ is a coefficient of the x/ term of 

d-3 

G(x) = 7 (x @ ath, 
{j= 60 


33. The method of claim 30 wherein said step of cal- 
culating said value E comprises the steps of: 
using said location L and a number a wherein 
0=a<d—1 to produce an index value for an ele- 65 
ment f(a,L) of said first table; 
retrieving said element f(a,L) of said first table refer- 
enced by said index value; 


wherein pre-computed constants Agy, Agc, Aad; 

Abe, Abd: Acd) Bab, Bac, Bad, Bac, Bad, Bea, Cad, 
Cac, Cad, Coc, Cod, and Ceq are functions of a, b, 
c, and d given by: 

Aab = (a4 @ a?) * (G'e_1 * G'd ® Ge * G'a_), 

Age = (074 @ a) * (G"5_1 * Gia ® G5 * G'g_t)s 

Aad = (a?4 @ a4) * (G"p_1 * G'e ® Gb * G"e—1), 

Abe = (a? @ a?) * (G"g_1 * G'g ® Ga * G'a_-1), 

Apd = (a2 @ 024) * (G"g_1 * Ge @ G"a * G"e_)s 

Aca = (0 @ a4) * (G"g_1 * G'h @ Gg * G'b—1)hs 

Bap = (a4 * a2 @ ab * 04) * (G"e_1 * G'a @ Ge * G'd-1)hs 

Bac = (a7 * a @ af * a9) * (G"p_1 * G"g © Gp * Ga), 

Bad = (a? * a4 @ at * a4) * (G"p_1 * Ge ® G5 * G'e—1)s 

Boe = (a? * a @ at * a2) * (G"y_1 * G'a © G's * Gat) 

Boa = (a> * 24 @ ad * a2b) * (G"g_1 * G'e @ G'a* G'em1), 

Bed = (a * 04 @ af * a7) * (G"g_1 * Gp @ G"a * G"b-1)s 


Cab = (a7 © a4) * (G"c_1 * Ga ® Ge * G" a1), 
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-continued 
Cac = (a? © a) * (G"b_1 * G'ag G G's * G"a_-1), 
Cad = (a9 @ a4) * (G"p_1 * Ge @ Gp * G'c_1), 
Coc = (a @ a) * (G"a_-1 * G"g @ Gg * G'a_1)s 
Cod = (ao @ a4) * (G"g_1 * G'e  G'a * G"e_1)s 
Cea = (0° @ a4) * (G"g_1 * Gb @ Ga * G"b_-1)s 
and G;’" is a coefficient of the x/ term of 
d—4 


Gx) = 7 (x Q a+) 


i= 


with G_1" defined as zero; and computing 
o1=N)/D and 02=N2/D. 


37. The method of claim 35 wherein t>4 and selected - 


coefficients Ra, Ry, Re, and Rgof said coefficients R;are 
each not equal to zero. 

38. The method of claim 35 wherein said step of cal- 
culating said values E; and E2 of said two data symbol 
errors comprises the steps of: 

computing non-zero parameters D, Nj, and No ac- 

cording to 


D = faliPfoL2) @ KOLiIYAGL2); 
Mv = Ra*fb,L2) @ Ro*fa,L2), 


Nz = Ro*fb.L1) @ Ro*haLy), 


wherein Ra, Ry, Re, and Rgare selected coefficients 
of said coefficients Rg and computing Ei; =N1/D 
and E2x=N2/D. 

39. The method of claim 35 wherein said step of vali- 
dating said locations and said values of said two data 
symbol errors comprises the steps of: 7 

calculating finite field logarithms of said values E, 

and E; 

testing a plurality t—2 of said remainder coefficients 

R; wherein i=4a, is4b, isc, and iséd, and further 

wherein each test comprises a sequential repeated 

block comprising the steps of: 

adding said finite field logarithms of said values E; 
and E2 to respective elements f(i,L;) and f(i,L2) 
of said first table to produce finite field loga- 
rithms of test values, 

calculating finite field antilogarithms of said finite 
field logarithms of said test values; 

calculating an EXCLUSIVE-OR sum of said finite 
field antilogarithms of said test values; 

comparing said EXCLUSIVE-OR sum of said 
finite field antilogarithms of said test values to 
said coefficient Rj 

correcting said two data symbol errors using said 
locations and values when each said EXCLU- 
SIVE-OR sum of said finite field antilogarithms 
of said test values is equal to said coefficient R;. 

40. The method of claim 28 wherein said step of ap- 
plying said remainder coefficients further comprises 
validating locations L; and values E; of a plurality e of 
symbol errors by testing a plurality of said remainder 
coefficients R; according to the equations 
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Ri = 2. o Ef fil). 
j=l 


41. The method of claim 28 wherein said step of ap- 
plying said remainder coefficients further comprises the 
steps of: 

computing a syndrome polynomial S(x) from said 

‘remainder polynomial R(x); 


10 : : 
generating an error locator polynomial o(x) from 
said syndrome polynomial S(x); 
locating errors using said error locator polynomial 
o(x); and 
15 evaluating errors using said error locator polynomial 


o(x) and said syndrome polynomial S(x). 

42. The method of claim 41 wherein said step of com- 
puting said syndrome polynomial comprises the steps 
of: 

(1) initializing a coefficient So and all other coeffici- 
ents S;of said syndrome polynomial S(x) to a coef- 
ficient Ro of said remainder polynomial R(x); 

(2) initializing a counter j to 1; 

(3) computing a finite field logarithmic partial result 
comprising a MODULO 2” — 1 sum of a finite field 
logarithm of a non-zero coefficient Rj and j*mo; 

(4) calculating a finite field antilogarithm of said par- 
tial result; 

(5) EXCLUSIVE-OR adding said finite field antilog- 
arithm of said partial result to said coefficient So; 

(6) MODULO 2”~—1 adding said counter j to said 
partial result; 

(7) calculating a finite field antilogarithm of said par- 
tial result; : 

(9) EXCLUSIVE-OR adding said finite field antilog- 
arithm of said MODULO 2”—1 sum and one of 
said coefficients S; 

(10) repeating said steps (6) through (9) for said coef- 
ficients S; wherein i=2 to MIN(d—2,11); 

(11) incrementing said counter j; and 

(12) repeating said steps (3) through (11) for said 
coefficients Rj wherein j=2 to d—2. 

43. The method of claim 41 wherein said step of gen- 
erating said error locator polynomial comprises the 
steps of: 

computing an nth discrepancy d, using sequential 
repeated blocks, each said block comprising calcu- 
lating a MODULO 2”—1 sum of said nth discrep- 
ancy d, and finite field product of a coefficient op; 
and a coefficient S,—; 

using a degree 1, of an nth error locator polynomial 
a(x) as an index into a table of software addresses 
of each said block; 

generating coefficients S;for i=11 to d—2 when said 
parameter is equal to twelve. 

44. The method of claim 41 wherein said step of lo- 
cating said errors when said error locator polynomial 
o(x) is of degree j greater than four comprises the steps 
of: ; 

evaluating said error locator polynomial o(x) for 
successive values of L until 
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65 
said step of evaluating comprising a sequence of 
repeated blocks, each said block comprising calcu- 
lating a finite field product of a~4£ and an EXCLU- 
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SIVE-OR sum of said parameter A and a coeffici- 

ent o7; of said error locator polynomial o(x); and 

maintaining a software address of a starting point for 
next said evaluation. 
45. The method of claim 41 wherein said step of lo- 
cating said errors when said error locator polynomial 
o-(x) is of degree j less than or equal to four comprises 
locating one of said errors using one of the steps of: 
when j is equal to four, calculating a finite field loga- 
rithm of a root of a quartic equation in a finite field; 
when j is equal to three, calculating a root and a finite 
field logarithm of said root of a cubic equation in a 
finite field; . 
when j is equal to two, calculating a root and a finite 
field logarithm of said root of a quadratic equation 
in a finite field; or ea 
when j is equal to one, calculating a finite field loga- 
rithm of a root of a linear equation in a finite field. 
46. The method of claim 41 wherein said step of eval- 
uating said errors comprises the steps of: : 
dividing said error locator polynomial o(x) by 

(x@a£) to produce a new error locator polynomial 

S(x) and calculating an error value E from said 

syndrome polynomial S(x) and said new error loca- 

tor polynomial (x), all in a single software loop 
comprising steps of: : 

(1) initializing a counter g=1, a remainder R=1, a 
denominator D=1, and a numerator N=o;7 

(2) calculating a MODULO 2”—1 sum of said 
remainder R and a finite field product of a finite 
field antilogarithm of said location L and said 
remainder R; 

(3) storing said MODULO 2”—1 sum as said re- 
mainder R and asa coefficient og of said error 
locator polynomial o-(x); 

(4) calculating a MODULO 2”—1 sum of said 
remainder R and a finite field product of a finite 
field antilogarithm of said location L and said 
denominator D; 

- (5) storing said MODULO 2”—1 sum as said de- 
nominator D; 

(6) calculating a MODULO 2”—1 sum of said 
numerator N and a finite field product of said 
remainder R and a coefficient S;.z of said syn- 
drome polynomial S(x); 

(7) storing said MODULO 2”—1 sum as said nu- 
merator N; 

(8) incrementing said counter g; and 

(9) repeating said steps (2) through (8) for values of 
said counter g up to and including j; 

terminating error correction unsuccessfully when R, 
D, or N is equal to zero after completion of said 
steps (1) through (9); 

calculating a finite field quotient of said numerator N 
and said denominator D; 

recording said a finite field logarithm of said finite 
field quotient as a parameter E’; 

calculating a finite field product of said finite field 
quotient and a finite field antilogarithm of —L*mo; 

recording said finite field product as said error value 

E; and 

adjusting coefficients of said syndrome polynomial 
S(x) using steps of: 

(a) initializing counter g=0; 


4,839,896 


5 


15 


20 


25 


45 


50 


55 


60 


65 


30 

(b) calculating a finite field antilogarithm of said 
parameter E’; 

(c) calculating a MODULO 27-1 sum of said 
finite field antilogarithm and a coefficient S; of 
said syndrome polynomial S(x); 

(d) storing said MODULO 2” —1 sum as said coef- 
ficient S; 

(e) calculating a MODULO 2”—1 sum of said 
parameter E’ and said location L; 

(f) storing said MODULO 2”—1 sum as said pa- 
rameter E’; ‘ 

(g) incrementing said counter g; and 

(h) repeating said steps (b) through (g) for values of 
said counter g up to and including j. 

47. In a decoder for an error detection and correction 
system using a Reed-Solomon code or related code of 
degree d—1 for detection and correction of a plurality 
of errors wherein a message block is comprised of N 
interleaved codewords of said code wherein codeword 
i is comprised of n;—(d—1) data symbols and d—1 
check symbols comprising a total of D data symbols 
stored in a data buffer means and N*(d— 1) check sym- 
bols stored in a remainder buffer means where d, i, N, 
and n;are positive integers, and further wherein the first 
check symbol in said remainder buffer means is remain- 
der coefficient Rg_2 of codeword D MOD N and the 
last symbol in said remainder buffer means is remainder 
coefficient Ro of codeword (D—1) MOD N, with other 
coefficients interleaved between, a method for access- 
ing said data buffer means and said remainder buffer 
means for détection and correction of said errors com- 
prising the steps of: 

(1) initializing a parameter DMN equal to said num- 
ber of data symbols D MODULO said number of 
interleaves N and initializing a counter I to zero; 

(2) if said parameter DMN is equal to zero, resetting 
said parameter DMN to said number of interleaves 
N; 

(3) computing a forward displacement within said 
remainder buffer means of coefficient Ro of code- 
word I by calculating N*(d—1)—DMN; 

(4) computing forward displacements within said 

" remainder buffer means of other coefficients R; of 
said codeword I by repeated subtraction of said 
number of interleaves N from said forward dis- 
placement of said coefficient Ro of said codeword 

ec 

(5) determining location(s) L; and value(s) E; of er- 
ror(s) in said codeword I; 

(6) computing a forward displacement within said 
data buffer means of a last data symbol within said 
codeword I as Fingx=d—DMN; 

(7) computing a forward displacement within said 
data buffer means of an error at a location L;>d—1 
as Fi=Fnax—N* Le 

(8) correcting said error at said forward displacement 
F; using error value Ej 

(9) repeating said steps (7) and (8) for all errors in said 
codeword I; ; 

(10) decrementing said parameter DMN and incre- 
menting said counter I; and 

(11) repeating said steps (2) through (10) for values of 


said counter I up to and including N—1. 
* * * * * 


